Sub-ballistic behavior in quantum systems with Levy noise 
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We investigate the quantum walk and the quantum kicked rotor in resonance subjected to noise 
with a Levy waiting time distribution. We find that both systems have a sub-ballistic wave function 
spreading as shown by a power-law tail of the standard deviation. 
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In the last decades the study of simple quantum sys- 
tems, such as the quantum kicked rotor (QKR) [l| and 
the quantum walk (QW) Q, have exposed unexpected 
behaviors that suggest new challenges both theoretical 
and experimental in the field of quantum information 
processing Q. The behavior of the QKR has two char- 
acteristic modalities: dynamical localization (DL) and 
ballistic spreading of the variance in resonance. These 
different behaviors depend on whether the period of the 
kick is a rational or irrational multiple of 4ir. For ratio- 
nal multiples the behavior of the system is resonant and 
the average energy grows ballistically and for irrational 
multiples the average energy of the system grows, for a 
short time, in a diffusive manner and afterwards DL ap- 
pears. Quantum resonance is a constructive interference 
phenomena and DL is a destructive one. The DL and the 
ballistic behavior have already been observed experimen- 
tally [HIB. On the other hand the concept of QW intro- 
duced in pi, 0] is a counterpart of the classical random 
walk. Its most striking property is its ability to spread 
over the line linearly in time, this means that the stan- 
dard deviation grows as o~(t) ~ t, while in the classical 
walk it grows as a(t) ~ t 1 / 2 . We have developed |,|| a 
parallelism between the behavior of the QKR and a gen- 
eralized form of the QW showing that these models have 
similar dynamics. In (Io| we have investigated the reso- 
nances of the QKR subjected to an excitation that fol- 
lows an aperiodic Fibonacci prescription; there we proved 
that the primary resonances retain their ballistic behav- 
ior while the secondary resonances show a sub-ballistic 
wave function spreading (<r(i) ~ t c with 0.5 < c < 1) 
like the QW with the same prescription for the coin [ll[ . 
Casati et al. 12j have studied the dynamics of the QKR 
kicked according to a Fibonacci sequence outside the res- 
onant regime, they found sub-diffusive behavior for small 
kicking strengths and a threshold above which the usual 
diffusion is recovered. More recently Schomerus and Lutz 
[lU investigated the QKR subjected to a Levy noise [l4| 
and they show that this decoherence never fully destroys 
the DL of the QKR but leads to a sub-diffusion regime 



for a short time before DL appears. 

In this article we investigate the QKR in resonant 
regime and the usual QW when both are subjected to 
decoherence with a Levy noise. In the case of the QKR 
the model has two strength parameters whose action al- 
ternate in a such way that the time interval between them 
follows a power law distribution. In the case of QW the 
model uses two evolution operators whose alternation fol- 
lows the same power law distribution. We show that this 
noise in the secondary resonances of the QKR and in 
the usual QW produces a change from ballistic to sub- 
ballistic behavior. This change of behavior is similar to 
that obtained for both systems when they are subjected 
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to an aperiodic Fibonacci excitation [1( 

Levy distribution - We consider two time step unitary 
operators Uq and U± in a large sequence to generate the 
dynamics of the quantum system. The time interval for 
the alternation of Uq and U± is generated by a waiting 
time distribution p(AT), where AT = iT with T a time 
step and i an integer. Then i is the number of times that 
Uq is repeated before U\ is applied once, e.g. the sequence 
of operators when the first interval is ATi = 4T and the 
second is AT 2 = 2T, is UxUqUqUiUqUqUqUq. In this 
paper we take p(AT) in accordance with the Levy distri- 
bution This distribution appears frequently in nonlinear, 
fractal, chaotic and turbulent phenomena [HI, [ItJ , it 
includes a parameter a, with < a < 2, and it is identi- 
cal to the Gaussian distribution when a — 2. When AT 
is large the asymptotic behavior of p(AT) is (l/AT) 1+a , 
this implies that the second moment of p(AT) is infinite 
when a < 2 and then there is no characteristic size for 
the temporal jump, except in the Gaussian case. It is 
just this absence of characteristic scale that makes Levy 
random walks scale-invariant fractals. As we are inter- 
ested in the asymptotic behavior of the QKR and QW, 
the most important characteristic of the Levy noise is the 
power law shape of the tail. To capture the essence of 
the Levy noise distribution, and simplify at the same time 
the numerical calculation, we define the waiting time dis- 
tribution as 
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Then, the mechanism to obtain the time interval AT in 
agreement with the previous discussion is the following: 



2 



a) we sort a stochastic variable 7 with uniform distribu- 
tion in [0, 1], b) we obtain £ as a solution of the equation 

7 = J p(t)dt and finally c) AT = iT where i is the integer 


part of £. 

Quantum kicked rotor.- The QKR Hamiltonian is 

p2 °° 
H= — + K cos 6^ 8(t-nT), (2) 
21 n=l 

where the external kicks occur at times t = nT with n 
integer and T the kick period, I is the moment of iner- 
tia of the rotor, P the angular momentum operator, K 
the strength parameter and 9 the angular position. In 
the angular momentum representation, P\l) = £H\i), the 
wave-vector is \^(t)) = J2"eL~ 00 a e(t)\£) an d the average 

energy is E(t) = = eEfc-oo^ 2 M*)| 2 > where 

e = h 2 /2I. Using the Schrodinger equation the quantum 
map is readily obtained from the Hamiltonian pj) 

00 

ai(t n+1 ) = Uija,j(t n ), (3) 

where the matrix element of the time step evolution op- 
erator J7(k) is 

U tj =i-U- £ >e- i W l J j - i (K), (4) 

J m is the 771th order cylindrical Bessel function and its 
argument is the dimcnsionlcss kick strength k = K/h. 
The resonance condition does not depend on k and takes 
place when the frequency of the driving force is commen- 
surable with the frequencies of the free rotor. Inspection 
of eq.(|4]) shows that the resonant values of the scale pa- 
rameter r = eT/2h are the set of the rational multiples 
of 47r, i.e. t = 47r p/q. In what follows we assume that 
the resonance condition is satisfied, therefore the evolu- 
tion operator depends on k, p and q. We call a resonance 
primary when p/q is an integer and secondary when it is 
not. 

With the aim to generate the dynamics of the sys- 
tem we consider two values of the strength parameter 
k\ and K2, and combine the corresponding time step 
operators Uq = U and Ui = U (k 2 ) in a large se- 
quence. We have proved in [l(| that the ballistic behav- 
ior is maintained in the primary resonances for any type 
of sequences of the operators Uq and Ui because the op- 
erators Uq and U\ commute. For the same reason the 
antiresonance p/q = 1/2 is not changed. Then we only 
need to study the secondary resonances in this work. 

Using the operators Uq, U\ and the waiting time dis- 
tribution ([1]) we obtain the wave function spreading as 
given by the exponent c of the asymptotic expression of 

the standard deviation a(t) = \Jj2iL-oo^ 2 l a ^)| 2 ~ 
The initial condition for the wave-vector is the position 
eigenstate |0), that is a (0) = 1. The average standard 
deviation o~(t) is numerically obtained running a code 



with an ensemble of 10000 trajectories for several thou- 
sands of T. The table © 

p/q 1/3 1/4 1/5 2/5 
c 0.87 0.92 0.59 0.89 ' 1 ' 

shows that c depends on the ratio p/q. It is calculated 
with Ki = 0.5, K2 — —0.5 and a = 1. The value of c re- 
mains unchanged when p/q is changed for (q — p) /q; this 
symmetry in c being a consequence of the trivial symme- 
try of the time step evolution operator U as was shown 
in [10(. We have verified that the exponent c has a de- 
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FIG. 1: The standard deviation for the QKR as a function 
of time in units of T, with the parameters ki = 1, «2 = — 1 
and p/q = 1/3. (a) a = 0.2 and c = 0.998; (b) a = 1 and 
c = 0.772; (c) a = 2 and c = 0.518. 

pendence with the strength parameters Kiand k 2 and its 
range is always between 0.5 and 1, thus the sub-ballistic 
behavior is maintained. In Fig.[T]a'(i) is plotted for sev- 
eral values of the parameter a, displaying the qualitative 
differences between the periodic case (a ~ 0, c ~ 1), 
the Levy noise case (0 < a < 2 and 0.5 < c < 1) and 
the Gaussian case (a ~ 2, c ~ 0.5). The exponent c is 
plotted in Fig. [2] as a function of a showing a clear depen- 
dence of c with a. We have also studied higher moments 
of order four and six, they all have smaller exponents 
than those obtained with a periodical sequence, thus the 
asymptotic behavior of these moments is consistent with 
the power-law behavior of the second moment. 

Quantum walk.- The standard QW corresponds to 
a one-dimensional evolution of a quantum system (the 
walker), in a direction which depends on an additional 
degree of freedom, the chirality, with two possible states: 
"left" \L) or "right" \R). The global Hilbert space of the 
system is the tensor product H s eg) H c where H s is the 
Hilbert space associated to the motion on the line and 
H c is the chirality Hilbert space. Let us call T_ (T+) the 
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FIG. 2: The exponent c of the power law of the standard 
deviation as a function of the parameter a. (a) The black 
dots correspond to the QKR with ki = 1, ki = — 1 and 
p/q = 1/3, the full line is an adjustment. The parameters for 
the QKR are.; (b) The white dots correspond to the QW with 
6*i = 7r/3, 82 = 7r/6, the dashed line is an adjustment. 

operators in H s that move the walker one site to the left 
(right), and \L)(L\ and \R)(R\ the chirality projector op- 
erators in H c . We consider the unitary transformations 

U(0) = {T-(g>\L)(L\ + T + <S> \R) (R\}o {I ® K{6)} , (6) 

where K(9) — a z e^' l9<7y , / is the identity operator in 
H s , and a y and a z are Pauli matrices acting in H c . 
The unitary operator U (9) evolves the state in one time 
step, + 1)) = U(9)\^(t)). Here we are generaliz- 
ing the QW to the case where different quantum coins 
are applied as in As for the QKR, we combine two 
different step operators Uq — U(6i) and U\ — U{9 2 ) 1 
with 9x ^ 9 2l into a large sequence where we apply the 
same Levy waiting time distribution. The wave func- 
tion spreading is given by the exponent c in a(t) — 

Jj2Zo i% (k(*)| 2 + IM*)| 2 ) ~ tc - We take as the ini- 
tial condition for the QW the position eigenstate |0), with 
chirality (1,0) for the calculations of Fig. [3] and Fig. [2J 
and (1, i)/v2 for the calculations of Fig. [H The standard 
deviation a(t) is numerically obtained for an ensemble of 
1000 trajectories for each value of the parameter a of 
the Levy distribution. In Fig. [3]cr(i) is plotted for differ- 
ent values of a and the sub-ballistic behavior is clearly 
shown. In Fig. [2] c is plotted as a function of the pa- 
rameter a, here a functional dependence between them 
is evident. In Fig.[2]c 1S plotted as a function of 9, where 
9 = 9\ = — 6>2This figure shows that the range of c is 
always [0.5, 1], thus the sub-ballistic behavior is indepen- 



FIG. 3: The standard deviation for the QW as a function of 
time in units of T with the parameters 9i — tv/3 and 62 = 7r/6. 
(a) a = 0.2 and c = 988; (b) a = 1 and c = 0.71; (c) a = 2 
and c = 0.546. 



dent of the value of 9 (except for the trivial cases 9 = 0, 
9 = 7r/2); additionally this figure is in concordance with 
Fig. 3 in ref. [TTJ] These figures show qualitative simi- 
larities with the corresponding figures for the QKR, and 
point to the parallelism between the QW and the QKR 
in the secondary resonance regime. Again in this system, 
the moments of order four and six are consistent with the 
power-law behavior of the second moment. Conclusion.- 
The quantum resonances and the DL of the QKR have 
been experimentally observed in samples of cold atoms 
interacting with a far-detuned standing wave of laser light 
[Til Il9l Ho] and in particular the secondary resonances 
have been recently observed by Kanem et al. [5j. On 
the other hand several systems have been proposed as 
candidates to implement the QW model. These propos- 
als include atoms trapped in optical lattices [2l| . cavity 
quantum electrodynamics 12211 and nuclear magnetic res- 
onance in solid substrates [23l [24| . All these proposed 
implementations face the obstacle of decoherence due to 
environmental noise and imperfections. Thus the study 
of the behavior of these systems, subjected to different 
types of noise, is very important for the design and con- 
struction of future technologies. Here we proposed the 
study of the QKR and QW subjected to noisy pulses with 
a Levy waiting time distribution. As Gaussian noise is a 
particular case of the Levy noise, then our study is open 
to wider experimental situations. We prove that for QKR 
and QW the Levy noise does not break completely the 
coherence in the dynamics but produces a sub-ballistic 
behavior in both systems, as an intermediate situation 
between the ballistic and the diffusive behavior. Then 
QKR and QW have essentially the same dynamical evo- 
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lution and our results fortify the previously established 
parallelism between them [U, . More generally we 

can say that the dynamical evolution of the QKR and 
QW show certain patterns that seem to be common to a 
greater class of systems that are defined mainly by their 
symmetries and not by their microscopic details. The 
existence of an universality in the behavior of these sys- 
tems suggests that one is allowed a larger flexibility in 
the choice of the physical systems to build quantum com- 
puters. It is important to highlight that the sub-ballistic 
behavior obtained in this paper for the QW and the QKR 
is essentially the same as that obtained for these systems 
when the perturbation follows a Fibonacci prescription 
[Tol . [IlT | . The reason may lie in the fact that behind the 
Fibonacci prescription hides the lack of a typical scale in 
the sequence of the operators Uq and U\ which leads to 
a power law [25j in the standard deviation, in the same 
way as for the Levy noise. 



FIG. 4: The exponent c that characterizes de power law of the 
standard deviation for the QW as a function of the parameter 
9. The evolution is obtained with 9 = 9i = —92 in the interval 
[0, 7r/2). The dots correspond to the calculation and the full 
line is a polynomial adjustment 
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